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We present the consistent SN formalism for curvature perturbations in anisotropic cosmological 
backgrounds. We employ our 5N formalism to calculate the power spectrum, the bispectrum and 
the trispectrum in models of anisotropic inflation with the background gauge fields in Bianchi I 
universe. Our results coincide exactly with the recent results obtained from in-in formalism. To 
satisfy the observational constraints the anisotropics generated on power spectrum are kept small but 
large orientation-dependent non-Gaussianities can be generated. We study the Suyama-Yamaguchi 
inequality for the amplitudes of the bispectrum and the trispectrum in the presence of anisotropic 
shapes. 



I. INTRODUCTION 

Recently there have been many interests in anisotropic inflation. This is partially motivated from the cosmological 
observations indicating some statistical anisotropics on cosmic microwave background (CMB) Jj. Although the 
statistical significance of the violation of statistical isotropy is not high 0, [sj , but nonetheless the possibility of having 
statistically anisotropic seed perturbations are intriguing. One useful method to quantify the statistical anisotropy 
is to vifrite the curvature perturbation power spectrum Vq in Fourier space for mode A: as Q 'P^ = "^0(1 + g*{k-n)^) 
in which h is the preferred direction in sky. Constraints from CMB and large scale structure indicate that < 0.4 
ii- 

The best method to generate anisotropic perturbations is to employ gauge fields during infiation. However, due to 
conformal invariance of C/(l) gauge fields in FRW background, the background gauge field energy density and their 
perturbations will be quickly diluted during inflation. Therefore, one has to break conformal invariance such that 
the gauge field energy density will not decay during inflation. One popular method is to consider a time-dependent 

gauge kinetic coupling such that the U{1) action has the form AC = —^^^^Ff^i,F^'^ in which 4> is the inflaton field 
and F^i/ is the U{1) gauge field strength. Furthermore, in order for the gauge field perturbations to be scale-invariant 
one has to choose /(</>) oc in which a{t) is the scale factor. These models in the context of anisotropic infiation 
and magneto-genesis were studied in great details in literature l7|-[42l[ . 

An interesting model of anisotropic inflation was proposed in [16| in which with f{(j>) oc the inflationary system 
admits an attractor solution in which the gauge field energy density, i.e. the electric field energy density, and the 
metric anisotropy reaches a small but cosmologically observable level. The cosmological perturbations for this model 
were studied in [43r}48j . Similarly, the cosmological perturbation analysis for a different model of anisotropic infiation 
(3l| involving a complex inflaton fleld charged under the U{1) gauge fleld were performed in I49|. These analysis are 
based on standard in-in formalism which proved technically difficult due to anisotropic background. On the other 
hand, experiences with SN formalism [50. - .57j in models of infiation with scalar fields showed that SN analysis are 
technically much easier to handle when calculating the curvature perturbations and their correlations such as power 
spectrum and bispectrum. This is specially the case for models in which there are not much interactions when the 
modes of interest leave the horizon and physically interesting effects, such as non-Gaussianities, originate from local 
type interactions once the modes are outside the horizon. Therefore it will be very helpful to extend the standard 5N 
formalism (50l - l56| to models of anisotropic backgrounds. This is one of our main goal in this work. 

There have been works in the literature employing the conventional 6N formalism for the models containing vector 
or gauge fields but the effects of anisotropic background were not taken into account, i.e. the gauge field is treated on 
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the same footing as the scalar fields in an FRW background. In this work we present a consistent 5N formalism for 
anisotropic backgrounds such as in [161] in which the background metric is in the form of Bianchi I. After presenting 
our 5N formalism we calculate the power spectrum and reproduce exactly the results in [i^, \^ . We also calculate 
the bispectrum which coincides exactly with the results of . 

Planck is expected to release its data soon. Any detection or otherwise of primordial non-Gaussianities from 
Planck will have significant implications for inflationary model buildings. Simple models of inflation predict almost 
scale-invariant and almost Gaussian perturbations. Therefore, any detection of primordial non-Gaussianity will go a 
long way to rule out or classify different inflationary scenarios. Non-Gaussianity may take different shapes in different 
models, for a review see [5^[59|. In models of inflation based on scalar fields the shapes of Bispectrum and Trispectrum 
are statistically isotropic. However, in models of anisotropic inflation, one obtains new shapes which are anisotropic. 

As an important consistency condition for single field inflation, a detection of local form bispectrum in the squeezed 
limit can rule out all single field models of inflation provided the system reaches the attractor solution [gO, [ll| so one 
can neglect the evolution of curvature per turbations on super-horizon scales and the curvature perturbations have 
the initial Bunch-Davies vacuum state j62l 163^ . As a different consistency condition, the Suyama-Yamaguchi (SY) 
inequality (63 |. [65l |. [66ll67j between the amplitude of the Bispectrum in the squeezed limit, Jnl, and the amplitude of 
the trispectrum in the collapsed limit, tnl, are expected to hold generally in models of inflation based on scalar fields. 
It is an interesting question to see if the SY inequality holds when the primordial perturbations are not statistically 
isotropic. We will study this question in the context of anisotropic inflation. 

The rest of the paper is organized as follows. In Section [ll] we present our 6N formalism. In Section IIIII we 
study a model of anisotropic inflation which provides a non-trivial set up to employ our SN formalism. In Section 
IIVI we present the bispectrum and the trispectrum analysis for the anisotropic inflation model and study the SY 
inequality. The conclusion and discussions are given in Section|Vl We relegates some technical details into Appendices. 

Note added: While this work was in its final stages the paper [68] appeared which has overlaps in Bispectrum 
and Trispectrum analysis in Section HVl with this work. 

II. 5N FORMALISM FOR ANISOTROPIC BACKGROUNDS 



In this section we extend the SN formalism [50|-|57[ to anisotropic backgrounds. First we present the background 
fields equations. After presenting the general metric perturbations, we look into the linearized fields equations and 
consider a systematic gradient expansion in terms of parameter e defined via 

e=— , (2.1) 

in which k represents the wave number in Fourier space. We demonstrate that the separate universe picture works, 
that is, in the limit e ^ 1 the background fields equations are locally hold inside each homogenized patch. 

A. Background Equations 

Our background is the the Bianchi I metric with the scale factors ai (t) , 02 {t) and 03 (t) 

ds'^ = -df + ai{tfdx^ -f a2{tfdy'^ + asitfdz^ . (2.2) 
We adopt the notations used in j69|] in which 

3 
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H.{t) = - , ff^i^iJ., (2.3) 



i=l 



in which is the Hubble expansion rate for the i-th spatial direction, i — 1,2,3 and a dot indicates the derivative 
with respect to t. 

The components of background Ricci tensor are 

R% = iii + Y^Hl (2.4) 

k 

R\ = (2.5) 

(5*j [h, + iHH^ . (2.6) 
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The background Ricci scalar is 

R^6H + 18H'^ - HkHk' . (2.7) 

k>k' 

To solve the Einstein fields equations we have to specify our energy momentum tensor. The general energy momentum 
tensor T^^ for an imperfect fluid has the form (70| 

Tf_,i, = {p + p)u^u^ +pgf^iy + q^iU^ + u^q^ + 7r^i/ (2.8) 

supplemented with the following conditions 

Qf^ = , Tt'^p = , TT^y = TT^p , TT^^ u" = , 

where u'^ is the fluid's four-vector velocity, p is the relativistic energy density, p is the isotropic pressure, tt^i, is the 
trace-free anisotropic pressure (stress) and q'' usually is referred to as "heat conduction", which is also the energy 
flux relative to u^. 

The special case of a perfect fluid is identified with tt^^ = = so we recover the standard form of r^,y for the 
perfect fluid 

Tf_cu = {p + p)u^,u^ +pgi^,i, (perfect fluid) . (2.9) 
For the comoving coordinate associated with the fluid we have 

w^ = (l,0) , u^-(~l,0), (2.10) 
so the Einstein equations can be read as 

m^^J2H,H, = ^ (2.11) 

f\ = q,^0 (2.12) 
M|iJ, = -3Mj.HH, + ^{p-p) + n\ (2.13) 

Here we have used the convention that Hi represents the background Hubble expansion rates while p,p and so on 
represent the background fluid's properties. We also defined H as the effective Hubble expansion rate appearing in 
Freidmann equation, Eq. (12. lip . Note, H here should not be confused with the Hubble expansion rate defined for 
conformal time usually used in literature. 

Finally, the energy conservation equation u^V ^T^'^ — results in 

- WpV^T^" ^p + 3H{p + p) + HjTt'jSf = 0. (2.14) 

in which again we have H = Hi/ 3. 

Note that, in this model, the Hubble parameter appearing in (0,0) component of Einstein equation, H, and the 
Hubble parameter appearing in continuity equation, H, are not equal. The difference between them is given by 

As a result H >%. 



1. Example: (7(1) gauge fields in an expanding background 

As an example of non-perfect fluid with anisotropic pressure and heat conduction, consider the standard U{\) gauge 
field theory in an expanding background. This theory will be the base of anisotropic inflation in next section. The 
action is Lem = —Ff^i^F'^" /A in which = — d„Afj is the field strength associated with the U{1) gauge field 

The electric field, i?^, and the magnetic field, i/^, are given by. 



(2.16) 
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and 



The electromagnetic energy-momentum tensor, T^Jj^, is 

rpp,v „ pail w V _ 
era a 

For an observer comoving with the fluid Tj^^^ can be written as [71j 



(2.17) 



(2.18) 



(2.19) 



where ri'^'^"^ is the four-dimensional totally antisymmetric volume element (770123 = \/— det g), h^i, = g^i, + Up^u^ is 
the projection matrix, = E^^E^ and = Hf^H'^, respectively, are the magnitudes of the electric and the magnetic 
fields and tt^^ is a traceless and space-like symmetric tensor given by 



<m = ^i^ + ^ ) ^ ^^^^ ^ Hf'H''. 



Eq. (j2.19p can be compared with the energy momentum tensor for a generic imperfect fluid defined in Eq. 
yields 



Pem 



+ H^) , 



6 



(2.20) 
which 

(2.21) 
(2.22) 
(2.23) 
(2.24) 



B. Perturbations 



Let us now consider the fields equations with perturbations. In our 5N analysis we adopt the notation used in 
[56j . The order of spatial derivative or the so-called gradient expansion is denoted hy e — k/ aH while the order of 
smallness of perturbations are denoted by 5. In principle, one has to consider different gradient expansion parameters 
ti for different directions ti = k/aiHi. However, to simplify the analysis we assume ~ e so there is no hierarchy for 
gradient expansions along different directions. 

The most general form of metric perturbations in this background is 



/ 



-2A 



\ 









2a?Vi 


1 9 


01037^^ 




2a|V'2 


aia27^^ 






2a§^3 / 



(2.25) 



in which A,P'^,ipi and 7*^ are small perturbations at the order 0{S) with 6-^1. Since in general there is no spatial 
symmetry left, we do not have vector or tensor perturbations. 

An important step in dealing with the gradient expansion ordering of Einstein equations is the order of the shift 
functions We note that at the background level /3* =0. As a result one expects that the background metric should 
be valid globally in the limit e — >■ and, as employed in (s^l, one can assume 



(2.26) 
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The ordering of /3* in Eq. (j2.26p was also obtained in [s^ with the assumption that the anisotropic pressure is first 
order in gradient expansion. We look into ordering of /3* more rigorously in Appendix [B] and verify Eq. (|2.26p . 

Now we have all the necessary materials for performing the gradient expansion analysis for perturbed Einstein 
equations. 

The (0, 0) component of perturbed Einstein tensor is 

(5G% = QAfL^ - H.ijj (2.27) 

in which H is defined in background Friedmann equation, Eq. ()2.1ip . Combining Eq. ()2.27p with the background 
(0,0) component equation, Eq. (|2.1ip . yields 

Ml, -6An'^ + Y,Hi'>p}j = p + Sp + 0{e'^) (2.28) 

Locally, as a function of (x, t), the above equation takes the form 

3M|,-H2(x, t) ^ p(x, t) + 0(e2), (2.29) 

in which 

H2(x,i) = ^H,(x,t)iJ,(x,i), (2.30) 

i>j 

with the following generalization of local Hubble expansion parameter Hi(x,t) 

H.(x.,)= ^-"',;y'" - (2.31) 

As a result one can readily associate the average local Hubble expansion rate i?(x, t) as 

H[x,t) ^ -2_^Hi{x,t) = ^^-j (2.32) 

i 

in which the background average Hubble expansion rate H is H = J^i Hi/3. 

Now we look at the energy conservation equation in its contracted form u^jV^T^^'^ — 0. At the background level 
the energy conservation equation is given by Eq. (|2.14p . Defining the fluid's proper time r via ^ = u^V^ ~ jip^ J^j 
the perturbed energy conservation equation is 



+ 3H(x, t) (p(x, t) + p(x, t)) + 



O(e'), (2.33) 



in which H{K,t) is the average local Hubble expansion rate defined in Eq. (|2.32l) . By using Eq. (IB2ip and (|B18p the 
above equation takes the following simple local form 

+ 3i/(x, t) ( p(x, t) + p(x, t))+y tt\{k, t)H,{K, t) = 0{e\,S^,Se) . (2.34) 

i 

in which 7r%(x, t) = 7f% + (57r%(x, t). 

So again we conclude that our separate universe recipe works and it is enough to replace any background function 
/(i) by its local form /(x, i) and also using new local directional Hubble parameters iJi(x, t). Our prescription will 
be satisfactory if we can also check the (i = j) components of Einstein equations which are identical to the dynamical 
equations of tt*. The diagonal spatial components of Ricci tensor to leading order can be read as 

R\ = (iJ, + 3HH^){1 - 2A) +4>, + 3Hij, + i^^ ^ i^j - H^A + 0{e^, 6^) . (2.35) 

j 

One can simply show that the above equation can be simplified as follows 

R'^ = + 3i7(x, i)iJ,(x, t) + 0{€\ 5^) , (2.36) 

dr 
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so the {i = j) components of Einstein equation simply modifies the corresponding background equation, Eq. (j2.13l) . 
as follows (for the off-diogonal components of Einstein equation see Appendix [C]) 

^^2^ dHjK,t) ^ ^SMlHH, + ^{p{^,t)-p{^,t))+7r\{^,t) (2.37) 

Now we have a complete set of local fields equations, Eq. (|2.29l) . Eq. (|2.34p and Eq. (|2.37l) . mimicking the 
corresponding background equations, Eq. (|2.1ip . Eq. (|2.14p and Eq. (|2.13p . with the local Hubble parameters 
ifi(x,t) defined in Eq. (|23T1) . 

The separate Universe approach discussion is now complete. The 6N formalism is also in hand noting that from 
thee above equation one has 

Ni{x,ti,t2) = / i?,(x,t)(l + = / Hidt+ / ip,dt (2.38) 

So one readily finds 

N,{^,tiM) ~ N,{t) = ^j,{t2) - Mti) (2.39) 
Now defining the average expansion by 

iV(x,ti,t2) - ^ ViV,(x,ti,t2) = f'H{^,t){l + A)dt (2.40) 

one obtains 

SN{^,t,,t2) = Af(x, t,,t2) ~ Nit) = V(i2) - i^ih) (2.41) 
in which i/'lx, t) is defined as the average of ip 

V;(x,<)^i^V^,(x,t). (2.42) 

i 

We are interested in curvature perturbation on surface of constant energy density. As was demonstrated in Appendix 
1X1 the curvature perturbation ^ defined in Eq. (|A9p via 

-C^ij~-Sp, (2.43) 
P 

is gauge invariant. The relation between C and SN therefore is 

a^,t)=SN{^,U,tf) (2.44) 
in which the initial surface is a flat surface -0 = and the final surface should be a uniform energy density surface 

Here a comment is in order. The diagonal components of the anisotropic pressure, Snl (no sum over i), are non-zero 
at the background level so their perturbations are expected to play some roles in the curvature perturbation analysis. 
However, the non-diagonal spatial components of anisotropic pressure and the heat conduction terms are absent at 
the background level so their perturbations will dilute quickly. The diagonal anisotropic pressure plays two different 
roles in the curvature perturbation analysis, a direct effect and an indirect effect. The direct effect can be seen from 
the continuity equation, Eq. (|2.34p . in which Stt} contributes to the Hubble expansion rate. This effect, by using 
Ea. (|2.34p . can be quantified as follows 

1 r*^,^ p(x,i) 1 



7V(x,i„i/)=/ H{^,t)dT^--J^ _ _ / d<Q(x,0, (2.45) 



in which Q{x.,t) is defined as 

Q(x,t) = 

p + p 



V — g^ + V^g^-u^.V^TT^'- 



7 



The above equation shows that the diagonal anisotropic pressure components ^7r^^ contributes to SN through their 
effect on continuity equation as captured by the term containing Q in Eq. (|2.45|) . 

The indirect effect of anisotropic pressure is more subtle and sometimes can be more important than the contri- 
bution from the term containing Q above. This effect can be understood as the back-reactions of fields responsible 
for anisotropic pressure on the dynamics of other background fields such as the inflaton field. The SN formalism 
automatically includes this indirect effect. We will see this effect in next section in application of our SN formalism 
for models of anisotropic inflation. 



III. ANISOTROPIC INFLATION 



In this section we present the model of anisotropic inflation with a U{1) gauge field originally presented in [4£ 
which provide a non-trivial setup to employ our SN formalism. 
The action is given by 



S = I d^x. 



M 



1 



(3.1) 



in which cj) is the inflaton field and F^i, — df^^Ai, — duA^ is the field strength associated with the U{1) gauge field A^. 

To employ the SN formalism, as usual we need to have a good control of the background dynamics. We assume 
that the gauge field has a non-zero classical value along the x-direction so A^ = (0, 0, 0). As a result, the 

background space-time is in the form of Bianchi I Universe with the metric 



= -df + a{tfdx^ + b{tf{dy^ + dz^) 



(3.2) 



In this view H = a is the average Hubble expansion rate, Ha = a/a and Hh = b/b are the expansion rates along the 
spatial directions x and y and &/H = {Hi, — Ha)/H is a measure of anisotropic expansion. 



A. The background dynamics 



The fields equations are given by 

(0)e"+4'^i,) = (3.3) 

+ 3a(/. + y^-/(0)/40)i2e-2a+4. ^ (3.4) 

1 ^2 ^ ^(^) ^ 1 ^2(^)^2^-2„+4. ^ 3^^2 (^2 _ ^2) (3 5) 

V{^) + ^-fiMle-^"^'" = Ml {a + 3d2) (3.6) 

i/2(0)i2e-2"+4- ^ Ml [Zaa + a) , (3.7) 

in which a dot indicates derivative with respect to t. 

The equation of motion for Ax (the Maxwell equation) is easily solved as 

^ = /(0)-'e-"(*)-4-(V, (3.8) 

where pA is a constant of integration. 

We are interested in the small anisotropy limit, \ct/H\ <C 1, so the background expansion is mainly supported by 
the isotropic potential term as in conventional models of inflation. In order for the anisotropy to be small, we demand 
that i? <C 1 in which 

In this view R measures the ratio of the electric field energy density, Pemi associated with the gauge field to the total 
potential energy density. Therefore, to have small anisotropics, we require pem ^ V . 
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Although the anisotropy is small, i? <C 1, so the Hubble expansion rate in modified Friedmann equation p.5p is 
mainly dominated by the isotropic potential term, but the back-reactions of the gauge field on the inflaton field induce 
an effective mass for the inflaton as given by the last term in Eq. (|3.4|) . This in turn will affect the dynamics of the 
inflaton field. As shown in 16] with the appropriate form of f{(j)) the the system reaches an attractor solution in 
which R reaches a subdominant but nearly constant value. For R to be constant, we need f{(f>) oc a" with n ~ —2. 
Indeed, the background expansion is given by 



a cx exp 







[-/ 


d<f>— 



(3.10) 



So if one chooses 



/ oc exp 



V 



(3.11) 



this yields / oc a". The exact form of / therefore depends on V{(f>). For the chaotic potential used in [la| we have 



V = -m^ 
2 



f{<P) = exp 



2Ml 



(3.12) 



with c a constant very close to unity. In our discussion below we take the form of /, in terms of a = e", to be 

in which and ??£ represent the value of the scale factor and the conformal time at the end of inflation. 
As shown in 16] the system reaches the attractor solution in which R is given by 

R = ^eff = \leH , (3.14) 



where we have defined / = and €h ^ H / H is the slow- roll parameter. Combined with the definition of R in Eq. 
we obtain 



'2a 



lenV. 



(3.15) 



As we shall see below, this equation will be the key equation to find 5N in terms of 5<j) and 5 A. 
Furthermore the anisotropy in expansion is given by 



(7 

a 



3 



During the attractor phase the inflaton evolution is given by 

P da- 



V 



(3.16) 



(3.17) 



Interestingly, this means that the back-reactions of the gauge field on the inflation field change the effective mass of 
the inflaton field as given by the second term above. 

Using Eq. p. 121) in Eq. p.l7p results in the following equation 



4M|,a(l - /) 



(3.18) 



in which 0e is the value of 4> at the end of inflation. We choose the convention such that ag = 0, so during inflation 
a < 0. Eq. p. 181) clearly shows the effect of the gauge field back-reactions on the evolution of the inflaton field. The 
fact that the evolution of the inflaton flcld is affected by the gauge fleld, as given by the correction factor (1 — /) in 
Eq. p.l8|) is the key to calculate 6N in the presence of gauge field. In passing, we comment that in the previous 
applications of 5N in the literature for models with the gauge fields, this important effect is not taken into account. In 
other words, 5N in these papers have been written with treating 5A^ in the same footing as 5(j) in an FRW background 
without taking into account the back-reactions of the gauge field in the evolution of inflaton fleld and in the dynamics 
of the anisotropic background. 
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In connection with our discussion in previous section the energy density, pressure, momentum density and stress 
associated with the electro-magnetic field are given by 



Vera 



(3.19) 

(3.20) 

(3.21) 
(3.22) 



and 





-2IeHH^ 

lenH"^ 










(3.23) 



Plugging the value of pem vato definition of R and using the attractor value Eq. p.l4p we obtains pe 
advertised before. Also Eq. p.l9p indicates that E = ^/sTchH and B = 0. 



RV as 



B. 5N in anisotropic inflation 

Our goal here is to calculate the curvature perturbations in this model by employing our 6N formalism. As we 
argued before the contribution of the gauge field into the Hubble expansion rate and total energy density is sub- 
dominant. This means that the surface of end of inflation is controlled only by the inflaton field. However, the gauge 
field plays an important role in Klein-Gordon equation and in the evolution of the inflaton field as can be seen in 
Ea. (l3T8l) . 

Perturbing Eg. (13. 181) we have. 



As a result 



= -ASN + ANSI . 



SN ^ ~^<j)S(l) + NSI . 



(3.24) 



(3.25) 



Note that, in order to connect to the standard notation we made the replacement a ^ N so Sa = SN from now on. 

The first term in Eq. (13.251) is the contribution of the inflaton field, while the second term is due to the back-reaction 
of the gauge field on the infiaton dynamics. Also by perturbing Eq. p.lSp we have 



€hSI = -12R6N + 4R 



(3.26) 



Now combining Eq. (I3.25|) and Eq. p.26p we have 



1 + ]6N = - 



2Ml 



ARN 5A,^ 
en Ax 



(3.27) 



Now using Eq. p.l4p we have RN/ch ^ IN. As we shall see, we require NI <C 1 in order not to produce too much 
anisotropy in power spectrum so we can neglect the second term in the left hand side of Eq. p.27p and 



6N : 



2Ml 



21 N 



SAx 



(3.28) 



This is our result for SN to linear order in terms of Scj) and SA. Interestingly, since in this model the leading 
contribution into the anisotropic power spectrum comes from the electric field instead of the magnetic field, we see 
that in SN only SA and not SA appears. This should be compared with the conventional models of SN involving 
scalar fields 0/ in which Sc/jj and not dtScjij appears. This is because for light scalar fields S(j)i become negligible on 
super-horizon scales once the attractor solution has been reached. 
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To calculate the power spectrum and the higher order correlations, we have to know the behavior of outside 
the horizon. For this purpose, we have to solve the mode function for SA^ with the initial Bunch-Davies vacuum deep 
inside the horizon. As shown in (48| the canonically normalized gauge field quantum fluctuations are given by 

'^^ = j:J(0f^^^-^'^(f^)j' (3-29) 

in which 

V = axCt)Vx{k) + al{-~i)V^{k) . (3.30) 

Here a\{ k) and a| ( fc ) represent the annihilation and the creation operators and for A = ± represents the 
circular polarization with the properties k.e±{k) — 0,k x e*±(fc) = ^ik e±{k) , ex{^k) = ex (A;)*, normalized via 
e\{k) . ey{k') = Sxy and 



ex^i{k)elj{k) = 6ij - hkj . (3.31) 



A 

The mode functions satisfy the evolution equation 

Vxiky + (^fc2 _ y^(fc) ^ , (3.32) 

where the prime denotes the derivative with respect to conformal time drj = dt/a(t). For / given in Eq. (I3.13P the 
normalized gauge field mode function is the same as that of a massless scalar field in dS space with 

Using this form of the wave function and the attractor solution Eq. p.lSp one can easily show that on super-horizon 
scales 

{k>aH) (3.34) 

A Y V2l^ 

In particular, we see that on super-horizon scale, SA^/Ax is a constant. 

Now we are in the position to calculate the total effect of the gauge field in curvature perturbation (. From Eq. 
p.45p and Eq. p.28p we see that there are two different terms that encode the contributions of the gauge field in 
C- Eq. (|3.28p encodes the indirect effects of the gauge field on 6N originating from its back-reaction on inflaton field 
dynamics. However, the direct contribution of the gauge field in SN is encoded in term Q in Eq. (|2.45l) . As we shall 
prove below, the contribution from the Q term in the curvature perturbation is negligible and the leading contribution 
of the gauge field in curvature perturbation is from its back-reaction effects in Eq. (|3.28p . 

Calculating Q from Eq. (|2.46|) yields (note that in this model is proportional to the product of the electric and 
magnetic fields and since in this model the magnetic field is zero therefore there is no correction from g'^) 



1 



(P + P) 
2H^ 



(H,- Ha) I en 



(p + p) 

I^HfH, (3.35) 



where we have used {p + p) ^ (jp = 2H^eH and Hi, — Ha = HI en from Eq. p.l6p . 
Perturbing Eq. p.35p . we have 



5Q = 2IHe„SI 

2RSN + AR- 

A 



2IH i -12RSN + AR^] . (3.36) 
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Now integrating Eq. (I3.36|) over t we have, 



6Qdt = / 2HI ( -12R6N 



ti 




21 N -12RSN + (3.37) 



To perform this integral, we have assumed hat H and R are nearly constant in the slow-roh approximation. Further- 
more, 6Ax/Ax is also nearly constant as can be seen from Eq. p.34p . 

Eq. p.37p indicates that the contribution of Q in 5N is at the order of INR ^ NP . However, as we shall see, in 
order not to produce too much anisotropy we require IN^ < 1 so NP ^ 1 and we can safely neglect the contribution 
of Q in 5N. 

In conclusion, the only contribution of the gauge field in curvature perturbations comes in 5N as given by the 
second term in Eq. p.28p . As a result we have 

Q = 5N 

We are interested in curvature perturbation power spectrum defined via 

(4^Cfe^) = (2^)353(fc, + fc2)Pc(fci) , Vdk)^^Pdk). (3.39) 

We decompose the power spectrum into the isotropic part, Vq, coming from the (50 contribution in Eq. (j3.38p and 
the anisotropic power spectrum, AP, coming from 5 A in Eq. p.38p 

Vc^Vo + AV. (3.40) 

As usual the isotropic power spectrum is given by 

To calculate the anisotropic power spectrum we note that (50 and SA are mutually uncorrelated so {S(j)5A)\^, — 0. As 
a results 



= ^ -^N^smH 

2tt^ enkf 

= 24:IN'^Vosm^e (3.42) 
in which the angle 9 is defined via cos9 — n.k. Now comparing this with the anisotropy factor defined via 

rdk)=ra[l+g*{k.nf) ■ (3.43) 

we obtained 

5* = ~2AIN^ (3.44) 

Very interestingly this is the result obtained in [Zsl . using the standard in-in formalism. The advantage of using 
SN formalism is that we only needed to use the background attractor solutions with the information about SA at 
the time of horizon crossing. This should be compared with the tedious analysis employed in (45l . l48l l49j using in-in 
formalism to calculate AV. Physically, one expects that the SN method to be applicable in this model. The reason 
is that all the dynamics between the infiaton field and the gauge field are in the form of local interactions and the 
dynamics of modes are trivial inside the horizon and at the time of horizon crossing. 

As mentioned in [H, in order to satisfy the observational bound jg*] < 0.3, and taking iV = 60 to solve the 
horizon and the fiatness problems, we require / < 10~^. As a result c is very close to unity and i? ^ 1 as advertised 
before. 

As emphasized in [4^ if one allows TV to be too large then the accumulative anisotropics produced from the IR 
modes can become too large. Therefore, the total number of e-foldings should not be too large in this model. 
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IV. BISPECTRUM AND TRISPECTRUM 



In this section we calculate the bispectrum and the trispectruni in the anisotropic inflation model studied in previous 
section using our SN method and compare the results with the corresponding results in [i^ and [6^ obtained from 
the standard in-in formalism. As we shell see the results for bispectrum and the trispectrum are in exact agreements. 

As usual, in order to calculate the bispectrum and the trispectrum in 5N formalism, we have to expand 5N to 
higher orders in perturbations. The expansion of SN to linear order is given in previous section in Eq. p. 381) . Here 
we generalize it to second order. To this purpose, we perturb the attractor solution Eq. p.l5|) to second order in 
fields perturbations 



T 



2/,. 



S<f)- 



f 

25N{ 1 



2SA^ 



f 

2SA, 



5A 
T 



f A 



(4.1) 



This formula gives a relation between 51, SN and different powers of Scfi and SA. On the other hand, Eq. p.25p from 
the perturbation of the evolution of (j){a) gives a relation between 51, 5N and S^i which is valid to all orders in SI and 
SN. Now plugging back Eq. (|4.ip into Eq. (|3.25p and keeping the leading corrections from / ^ 1 we obtain 



SN = 

in which to leading order in / 



N^SA 



Na 



+ 



N 



N^j^S^SA, 



Na 



and 



N 



2Ml 



2IN 



A 



N^ 



N 



p 



2/. 



AfJ, ^ Ml f 



AA 



2IN 



N 



4IN 



<pA 



A f 



(4.2) 



(4.3) 



(4.4) 



Having calculated SN to second order in Eq. (|4.2p . we can calculate the bispectrum Bi^{ki, k2, k^) defined via 

(C(fci)C(fc2)C(fc3)) = {2TTfS''{ki+k2 + k3)Bc{h,k2,k3) (4.5) 

There are three contributions into bispectrum; (a): N^^^N^^N^^{S(j)*), (b): N ^^N ^N ^{SA"^) and (c): 
N ^j^N^^N j^{S(t>'^SA^). The term (a) is the one expected from scalar field theory and is very small. The term 

(b) is purely from the gauge field while the term (c) is from the mixing of inflaton and the gauge field. One expects 
that the contribution of term (c) to be sub- leading as compared to the contribution of the term (b). Indeed, a direct 
analysis shows that the ratio of (b) to (c) is N so for iV ~ 60 one can safely neglect the contribution from the term 

(c) . In conclusion, the leading contribution to the bispectrum comes from {SA'^) and 



(C(fcl)C(fc2)C(fc3)) 



]^N j,j^{k^)N j,{k2)N j,{k:i) J ^(M,(fci)M,(fc2)M.(p)M,(fc3 -P)) +2perm. 



4I^N{ki)N{k2)N{k3) 



(2^ 



{SA,{ki)SA,{k2)SA,{p)SMk:,-'p)) + 2pcrr 



(4.6) 



in which N(ki) represents the time when the mode ki leaves the horizon. Now in the Coulomb gauge Aq — 0, the 
gauge field perturbations 5Ai{k) are given by [il] 



SA, 



E 



V3H 



(4.7) 



Plugging these in Eq. (14. 6p yields 



(C(fci)C(fc2)C(fc3)) - 288IN{ki)N{k2)N{k3) C{ki, k2)Po{ki)Po{k2) + 2pcrm 



{27rfS^(^ 



fci + fc2 + fca , (4.8) 
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in which the momentum shape function C{ki, k2) is defined via 

C(fci,fc2) = - (Mi-nf - (k2.nf + (ki.n) (k2.n) (kiM)^ (4.9) 
To obtain Eq. (|4.8I) we have used Po(fci) — ik'^^^M'^ ^'^^ isotropic power spectrum and 

^ 2Ienk^Ml {^^^ ' ^^'^^^ j ^^^^ ^ K^' + ^^-^^^ 

Using Eq. (|4.8p . one can calculate the bispectrum as 

B^{h,k2,h) = 288/iV(fci)iV(fc2)iV(fc3) (c(fci,fc2)Po(fci)-Po(fc2) + 2perm.) (4.11) 

This completes our results for the bispectrum. As expected, the shape of the bispectrum is anisotropic. Very 
interestingly our formula Eq. (|4.8p and Eq. (|4.1ip agree exactly with the result of (48|J obtained using the standard 
in-in formalism. 

To calculate /nl we go to the squeezed limit fci <C fc2 — ^3 in which 

fNL{ki,k2,ks) = lim iM^ll^. (4.12) 

■' ^ ' ' ^' fci->ol2 PQ{ki)P(k2) 

In this limit we get 

/jVL - 240/7V(fci)7V(fc2)'C(fci,fc2) (fci<A:2~fc3) (4.13) 

~ 10iV|5,|C(fci,fc2) (4.14) 
Taking ^ 60 and \g.^, \ ^ 0.1 and neglecting the orientation-dependence in /jvl this leads to large non-Gaussianity 

fNL - 60. 

Now we are in the position to calculate the trispectrum of our model. The trispectrum is defined via 

(C(fci)C(fc2)C(fc3)C(fc4)) = {2nfS^ (fci +k2 + k3 + h) rc(fci, fc2, fca, h) ■ (4.15) 

In the collapsed limit ki + k^ — k2 + k/^ — we can calculate the parameter tjvl via 

rj^Kih)- lim 1 T,ik,,k2,hM 

^ ' k,+k,^o'^Pd\ki+k3\)Pciki)Pcik3) ^ > 

The trispectrum (C(fci)C(^2)C(fc3)C(^4)) have 6 contributions in the forms of 

and (SA^). As in the case of bispectrum, one can easily check that the last term has the dominant contribution 

in trispectrum. Therefore to leading order we have 

(C(fcl)C(fc2)C(fc3)C(fc4)> ^ N j,^{k^)N j,^{k2)Njks)Njk^) J ^ J ^ 

'6A.dk:i)5Adki)5A,{q)5A,{ki - q)SAj{p)SAj{k2-p)) + Sperm. 



3456/A^(A:i)A(fc2)iV(fc3)A(A:4) (Dik^, k^, ki + k:i)P{k^)P{ki)P{\ki + ki\) 
+ llperm. ) (2Tif 5^ (ki +k2 + h + h) , (4.17) 



in which D[kz, fc4, fci + ^3) refers to the trispectrum's shape function and is given by 

D{k^, k4, ki + ks) — 1 — (/c4.n)^ — (ks.n)^ — (fci + k^.Ti f + {k^.n)(k4^.fi){k^.ki) + (fc4.Ti)(fci + k^.n){ki + fc3.fc4) 

+ (fc3.n)(fcTTfc3.n)(fcT+~fc3-^3) - (fc3-^i)(^4-n)(fcT^ (4.18) 
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Comparing Eq. (|4.18l) with the definition of trispectrum we obtain 

Tc(ki,k2,k3,h) ^ M56IN{ki)N{k2)N{k3)N{h)( D{k^ ] .(4.19) 



Now going to the collapsed limit fci + fcs = fc2 + fc4 = and using the definition of tn l given in Eq. (|4.16p we obtain 

TNL{kuk2, fcs, h) ~ 3A56IN{k3)^N{k4)^D{k3,k4, ki + fcs) . (4.20) 

As in the case of bispectrum the trispectrum is anisotropic so ttvl has direction-dependence. Comparing our 
trispcetrum with the results of obtained from in-in formalism, we find the exact agreement between these two 
results. 

Now we are in the position to check the SY inequality between fj^^ and tjvl which states 

TNL > (^/a^l) • (4.21) 

The importance of SY inequality as a tool to rule out inflationary scenarios was studied in [s^ [6^ . The SY inequality 
presented in the form of Eq. (|4.2ip is for the models in which /^vl and t^vl are either scale-invariant or have the 
same scale-dependence. In our case, see also [zl, Izll) we have complicated shape-dependent for Jnl and tnl so the 
original SY inequality as given in Eq. (14.211) is not applicable. Instead, in ^3] a general integral representation of SY 
is proved which states 

J d\^d?q2TML[qi,k ~ qi,q2,-q2 - k)Pc{qi)PM > l^j dPq^NLiq.-q- k,k)P^iq)^ , (4.22) 

in which fc — > 0. As demonstrated in [gl] this integral from of SY inequality does hold in our model if we assume 
(7* < 1. Indeed, the condition < 1 is necessary from the observational constraints on the power spectrum and the 
consistency of our starting assumption that the anisotropic power spectrum is sub- leading so AV < Vq. 

To see qualitatively that the SY inequality in its simple form, Eq. (|4.2ip . does hold in our model, let us neglect the 
direction-dependence in /jvl and tjvl coming from D(k3, k^, ki + k^) and C'{ki, ^2). As a result 

mL{kuk2,k3,k,) ^ 1 ^^^^^^ ^^23) 



(|)2/jVL(fc3)/jVL(fc4) 9* 

in which we have taken N{ki) = N and ~ —2AIN^ from Eq. p.44p . Demanding that < 1 from the cosmological 
observations and also from the consistency of our analysis we conclude that tnl > ^{SnlY so the SY inequality does 
hold. 



V. CONCLUSION AND DISCUSSIONS 



In this work we have presented the consistent 5N formalism in anisotropic backgrounds. We have demonstrated 
that the separate universe approach works. In each homogenized patch the local continuity equation and the local 
Friedmann equation hold which have the same form as the corresponding background equations. We note that 
the Hubble expansion rate appearing in continuity equation, _ff(x, t), and the Hubble expansion rate appearing in 
Friedmann equation, 'H(x, t), are different. 

The anisotropic pressure has two different effects in 5N analysis. One is the direct effect encoded by the term 
containing Q in Eq. (|2.46p . The second effect is indirect and comes from the back- reactions of the source of anisotropic 
pressure on the dynamics of other fields, such as the inflaton field. We have calculated these two effects in model of 
anisotropic inflation containing a U{\) gauge field. We have shown that the second effect, i.e. the back-reaction effect, 
is much larger than the direct effect coming from the Q term in 5N formula. In previous works in the literature, this 
back-reaction effect during inflation was not taken into account. The gauge field contribution 5Ai is added trivially 
as a non-interacting field during infiation. 

Taking into account the back-reaction of gauge field on infiaton dynamics, we have calculated 5N to linear and to 
second order in perturbations. We have demonstrated that our 5N formalism exactly reproduces the power spectrum 
and the bispectrum results obtained in previous works using standard in-in formalism. This is a non-trivial verification 
of the validity of our 6N analysis. The advantage in using 5N formalism is that all we need to know to calculate the 
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power spectrum and higher order correlations is the background dynamics and the profile of gauge field fluctuations 
on super-horizon scales. This method seems to be considerably simpler than the standard in-in formalism. 

We also calculated the bispectrum and the trispectrum in anisotropic inflation model. The bispectrum and the 
trispectrum are both orientation-dependent and scale-dependent. As a result the SY inequality in its simple form is 
not applicable. However, a generalization of the SY inequality in its integrated form indeed holds. 

Large amount of non-Gaussianity with /nl ^ 60 can easily be generated in this model. It is an interesting exercise 
to compare the predictions of this model on CMB and large scale structure formation along [t^ . [75| . 
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Appendix A: Metric perturbations and Gauge invariant perturbations 

In this Appendix we specify the properties of metric transformations under the general coordinate transformation 
and construct the gauge invariant curvature perturbations (. 

The background metric is given in Eq. (|2.2I) . The most general form of the scalar perturbations for the Bianchi I 
metric is introduced in Eq. (j2.25p . For the later convenience we introduce new variables /3* and 7*^ as 

dj^ = p\ d.djf' = (Ai) 

with no sum on repeated indices. 

Consider the general coordinate transformation 

^ .^M + ^ ^ (^^0 ^ Q^^\^ (^2) 

in which and — diC for i — 1,2,3 are scalars. Under the coordinate transformation Eq. ()A2|) we have 

6gf,,y 5g^^ - g^,^^^ C - 9au d^,^ - 9,.^" (A3) 
in which is the background Bianchi metric given in Eq. (j2.2p . More explicitly, one can check that 

(A4) 
(A5) 

(A6) 
(A7) 

(A8) 

If we apply the gradient expansion approximation df = 0{e^), then Q defined via 

^ _ (^1 + V'2 + V^3) {Hi +H2+ H3) Sp _ ^Sp ^^g^ 
3 3/9/3 

is gauge invariant and can be interpreted as the average curvature perturbations in our setup. 

Appendix B: Gradient Expansion Ordering of Perturbations 

In this section we estimate the ordering of /?* and 5q^ and calculate the contributions of Sq^ and Sirf^ in the energy 
conservation equation, Eq. (I2.33|) . 



A 












at 










2dfe 




- 


_ _ 


0.1 
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First of all let us check the transverse conditions on the heat flow and anisotropic pressure. By definition one has 

M^q^ = , w^tt"; = 0. (Bl) 

Noting that /3 — 0{5) in which 5 as usual is the order of perturbations, the fluid's 4- velocity to the leading order can 
be read as 



1 



1 + A 







-I -A, 



1 + ^ 



(B2) 



From the background equations we conclude that q'' and tt^ are zero. Now using the transverse condition (jBip one 
concludes that to all order 



6qo = , 6n''o = 
This also yields Sq'^ — 0{S^). For the ordering of Sn^i one has 

Sir", = afSn'o + Ptiih = ^tt^ 



(B3) 



(B4) 



We will use this equation later in order to find the ordering of the gradient expansion of perturbations. 

Let us now look at the gradient expansion ordering of /?'. For this we look at 5G'^i and (5G'o components of Einstein 
equations. With some efforts one can show that 



(B5) 



(B6) 



The easiest way to see this is to adopt the local inertial frame in which F^^ = 0. Therefore, the corresponding Einstein 
equations yield 



ST\ = eO(6)+pO{S) 

5T\ = P'{R\~R''o)+tO{6)+pO{6). 



(B7) 
(B8) 



Similarly, using the spatial components of energy momentum conservation one can put limits on ST'^i. The continuity 
equation V^^T/^ to leading order yields 



{dQ^'iH-Hi)5T\-a^^H,5T\ = -H,l3, {T\, - f%) + eO{S) + f30{S) 

= -H,I3, {R\ - R\) + eO{5) + pO{5) . 

One can show that the expression, T*; — T^o = R^i — R^o is non-vanishing in general and is of the order of H . 
Plugging Eqs. (|B7|) and (|B8|) into continuity equation Eq. (jB9|) yields 



{B9) 



[do + iH - Hi)6T'', =eO{5) + pO{5). 



(BIO) 



This equation shows that ST'^i has decaying solutions approximately like 1/a^. So one can readily deduce that 5T'^i 
should be higher order in gradient expansion as 

5T°, ^e^O{6)+el30{6). (Bll) 

Before discussing about the consequences of the above equation it is more convenient to rephrase Eq. (|B8|) as follows 

= p^{f\-f\)+eO{5)+l30{5), (B12) 

Furthermore for (JT'^i we have 

5T\ = [p + p)fi,-5q, + 5^\ + eO{5)+pO[5) (B13) 

By using Eq. (jB12[) to eliminate Sq^ in favour of /?* and Eq. (|B4[) to express Sir'^i as a function of one can show 
that the leading order terms of /3' cancel each other and one obtains 



(B14) 
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On the other hand, comparing Eq. (jBlip with Eq. (jB14p . one obtains the following result for the ordering of /3* 

=0(e), (B15) 

This also yields 



5q^ = 0[e) 
5t:\ = Oie) 



(B16) 
(B17) 



Now we investigate the contribution of heat flow in continuity equation Eq. (|2.33p . Using Eq. (jB16l) and Eq. (|B17p 
one finds that 5q^ ^ and by noting that the background value of is also zero, we get 



(B18) 



As a result one can deduce that heat conduction can be ignored in the continuity equation at the first order of 
perturbations and gradient expansion. 

Now it is time to calculate the contribution of anisotropic pressure on the continuity equation 

w'^V^TT^ = u'^V,^% + ,5K9,7r% + u'^r[;^7r% - u'^r^,7r%) (B19) 

By noting that 5u^^ — [^^j 0] one has 

w^V,7r% = u^V,^% + AH,tt\ + {do + 3H)Sn% + d^ir'o - H,Stt\ ~ ^,Tr\ (B20) 



Noting that Sttq = this yields 



(B21) 



Appendix C: (i 7^ j) Components of Einstein equation 

In this Appendix we look into off-diagonal components of spatial Einstein equations MpSG^j — 6Tj for i j- These 
equations are trivial at the background level. The leading order perturbation equations lead 



7ij + 3i/7i. 



1^J 



-Stt, 



(CI) 



in which e'-''' is Levi-Civita symbol. As one can see the above equation has decaying solutions. This is due to the fact 
that the background metric does not admit off-diagonal spatial components. Besides that, note that to linear order 
in perturbation, these non-diagonal terms do not show up in other components of Einstein equations so "fij evolves 
independent of other metric and matter perturbations. Indeed, one can easily see that 7.^ can be eliminated by an 
appropriate choice of gauge, so the (i ^ j) components of Einstein equation do not carry information for our SN 
formalism and can be ignored consistently. 
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